Universal amplitudes of the Casimir-like interactions between four 
types of rods in fluid membranes 
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Abstract. - The fluctuation-induced, Casimir-like interaction between two parallel rods of length 
L adsorbed on a fluid membrane is calculated analytically at short separations d <^ L. The rods 
are modeled as constraints imposed on the membrane curvature along a straight line. This allows 
to define four types of rods, according to whether the membrane can twist along the rod and/or 
curve across it. For stiff constraints, all the interaction potentials between the different types of 
rods are attractive and proportional to L/d. Two of the four types of rods are then equivalent, 
which yields six universal Casimir amplitudes. Repulsion can occur between different rods for soft 
constraints. Numerical results obtained for all ranges of d/L show that the attraction potential 
reaches ksT for d/L ~ 0.2. At separations smaller than dc ~ L{L/lpY^^ , where Ip is the rod 
persistence length, two rods with fixed ends will bend toward each other and finally come into 
contact because of the Casimir interaction. 
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When a field with long-range correlations is confined 
between objects that constrain its fluctuations, generic 
long-range forces arise. The most famous example is the 
Casimir interaction between two uncharged metal plates, 
which originates from the constraints they impose on the 
quantum electrodynamical vacuum fluctuations [1]. Simi- 
lar interactions arise between objects immersed in a clas- 
sical fluid, due to reduced thermal fluctuations. Fluids 
with long-range correlations produce long-range Casimir- 
like forces with universal characteristics. These interac- 
tions have been predicted to occur in binary mixtures 
close to a critical point, superfluids, liquid crystals and 
fluid membranes or interfaces [UIJ]. Direct measurements 
of thermal Casimir-like forces have been achieved recently 
between a colloid and a surface in a critical binary mix- 
ture [1]. 

In biomimetic membranes, which are fluid interfaces 
characterized by their bending rigidity [5], Casimir- 
like forces have been predicted to arise between phase- 
separated domains [6], integral proteins [6], adsorbed stiff 
rods [7] and semi-flexible polymers [8]. All these inclu- 
sions impose constraints on the thermal fluctuations of 
the membrane shape. The interactions between inclu- 
sions within biomcmbrancs, and their organization, plays 



a central role in many cellular regulatory processes and 
is of fundamental importance [9]. Biologically relevant 
stiff polymers adsorbed on membranes include actin fil- 
aments [iniin], FtsZ rings [E] and the /3- Amyloid fib- 
rils involved in Alzheimer's disease [T31II1]. Adsorption 
of goethite (a-FeOOH) nanorods on membranes has also 
been reported [15j . 

In this paper, we study the fluctuation-induced Casimir- 
like interaction free energy J^c between two parallel rods of 
length L adsorbed on a mcmbrancQ, at short and interme- 
diate separations d < L. It can be calculated by studying 
rods whose only effect is to restrict the membrane fluctu- 
ations: then, J^c corresponds to the d-dependent part of 
the free energy of the system. The Casimir-like force can 
be obtained by direct differentiation of Tc- We focus on 
membranes dominated by their bending rigidity k, i.e., we 
assume that the length scale (k/ct)^/^, which compares its 
effect to that of the membrane tension a, is much larger 
than all other relevant lengths. Other cases are briefly dis- 
cussed at the end of this paper. In the asymptotic limit 
L, previous analytical calculations showed that Tc is 
attractive and decays as k'B,T{L/d)'^ [7], where fee is Botz- 
mann's constant, and T is the temperature. Thus, even for 



^Our model is valid both for embedded and for adsorbed rods. 



p-1 



A.-F. Bitbol et al. 



d ~ 3L, its magnitude sa lO^'^fceT is far too small to ef- 
fectively bring rods together. At short separations d <^ L, 
one expects a constant scale-free Casimir interaction per 
unit length, i.e., J-q oc —k^TL/d. Effective attraction at 
d < L will thus result if the numerical prefactor is of or- 
der unity. Golestanian |H| showed that the finite bending 
rigidity of semiflexible polymers, or polymer rods, embed- 
ded in membranes, has dramatic effects on their Casimir 
interaction at short separations: shape fluctuations, which 
locally reduce the separation of the polymer rods, become 
unstable. Precisely, the free energy eigenvalues of these 
modes were shown to be X{q) ~ {Kr/kBT)q' - {2A/-K)/d^, 
where is the bending rigidity of the polymer rods 
and A ~ 2.91. Instabilities occur if A((jniin) < with 
ftnin ~ TJ'/ij i-c, for d < dc, where 

Here, Ip = Hr/{k^T) is the persistence length of the rod. 
For actin filaments [lOlfTT] of length L ~ 100 nm, with 
£p ~ 17/^m, this yields dc/L w 5x IQ-^. For L/lp ~ lO"", 
which may be realistic for actin bundles, the undulation 
instability threshold goes down to dd L k, 10^^. It is thus 
realistic to consider stiff rods at separations in the range 
0.01 < d/L < 1. In this paper we shall give analytical 
results for d -C L, and numerical results for all ranges 
of d/L. 

Let us describe the shape of the membrane by its height 
z ~ h{r) above a reference plane parametrized by r = 
{x,y). We shall refer to the latter coordinates as "in- 
plane" and to z as "out-of-plane" . In the Gaussian limit 
of small deformations around the fiat shape, the effective 
Hamiltonian of the membrane reads [S]: 

n^jdr^iy^hf. (2) 

Note that the Gaussian curvature does not contribute, due 
to the Gauss-Bonnet theorem [TB] , and that the membrane 
tension a has been neglected as previously discussed. The 
presence of a rigid rod lying within the membrane (or ab- 
sorbed on it) is traditionally described by imposing the 
constraint |6l[7]: 

h{r)^a + b-r (3) 

for r in a narrow rectangle [xq — e, + e] x [0, L], where 
a and b describe translational and tilt degrees of freedom, 
respectively. This is fully equivalent to imposing in the 
same domain the constraints: 

hyy = , ft.^^ = , hxy ^0, (4) 

where hij = didjh. An infinitely narrow rod can be 
described by imposing conditions ([4]) for x = xo and 
y G [0, L]. The single condition hyy = 0, however, is 
sufficient to specify a rod. In the absence of the other 
conditions, the membrane can curve across the rod and 
twist along it. It follows that we may consider four types 
of rods (fig. [TJ: 



1. Curved-twisted (ct): hyy = 0, 

2. Curvcd-nontwisted (ct): hyy = h^y = 0, 

3. Noncurved-twistcd (ct): hyy ~ hxx ~ 0, 

4. Noncurved-nontwisted (ct): hyy = h^x = hxy = 0. 

The actual constraints imposed on the coarse-grained field 
/i by a given rod depend on its microscopic properties, 
such as its width and its interactions with the lipids that 
constitute the membrane. Let us introduce the natural 
short length cutoff a of the theory, which is of order of 
the membrane thickness (it cannot be smaller since the 
membrane is considered as a thin surface): the field h is 
constructed through a coarse-graining, for instance a sim- 
ple moving average, of the height of each lipid molecule 
over the lengthscale a. First, let us consider a rod with 
width smaller than a that interacts with the lipids through 
short-range (e.g. hydrophobic) interactions, which can be 
neglected at a distance a. Such a rod cannot affect h^x or 
hxy, and it will thus be of the "ct" type. On the contrary, 
a rod with a width equal to, or a little bit larger than a 
will be of the "ct" type. Indeed, this width is sufficient 
for the thin fiat region arising from the presence of the 
rod to remain after the coarse-graining. A rod which has 
strong anisotropic attractive interactions with the lipids, 
with a range at least of order a, would behave as a "ct" 
rod too. In the case of intermediate ranges of interactions, 
the rod could constrain the slope but not the curvature of 
the membrane, in which case one would obtain the "ct" 
type. Finally, the "ct" type could correspond to a rod 
with width of order a and with short-range interactions, 
but composed of units free to twist relative to one an- 
other. In our coarse-grained description, all such rods can 
be considered as one-dimensional since their width is of 
order a or smaller. Hence, the microscopic properties of 
a rod (width, interactions) will appear only indirectly in 
our work, through the rod type. 

Let us first calculate the Casimir interaction between 
two parallel rods of the simplest "ct" type, in the d ^ L 
regime. Let us assume periodic boundary conditions with 
period L along the y direction. Generalizing the constraint 
hyy = to rods with a finite rigidity, we account for the 
presence of the rods by adding to % the following interac- 
tion energies between each rod and the membrane: 

for i £ {1, 2}, with xi = and X2 = d > 0. Thus, we allow 
out-of-plane fiuctuations of the rods with a bending rigid- 
ity K± ~ k-QTl^. The hard constraint hyy = is recovered 
in the limit kj_ — > oo. Note that we implicitly assume in 
any case an infinite in-plane rigidity «;|| oo, which is jus- 
tified if d/L is not too small, as discussed above. Adapt- 
ing the field-theoretical method of Li and Kardar [17] to 
these general constraints, we introduce two auxiliary fields 
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Casimir interactions between various types of rods 




where 



Fig. 1: Four types of rods for which the membrane may or may 
not curve (c) across the rod or twist (t) along it. 



ipi{y), and we make the substitution 



fdyi^^ip^+i f dr h{r)rp'/ {y)6{x—Xi) 



(6) 



in the partition function Z ~ J Vh exp[— /3('H + Ti'^ + H2)] 
of the membrane containing the two rods, where /? = 
{k^T)~^. Performing the Gaussian integral over h{r) 
yields 



(7) 



where ip^{q) = Vi(2/) cxp(-igy), = (-01(9) i'2{q)Y 



and 



M{d,q) 



r(d,g) ei^ + r{o,q) 



(8) 



kBT 4 dp 
kBT 



|g|(l + |g|x,)e-l«l-' 



(9) 



is the Fourier transform along y of T{xi,y) = 
d*G{xi,y)/dy^, where G(r-r') = {r\{l3K\7^)-^\r') is the 
correlation function of the free membrane. Performing 
the Gaussian integrals over the fields tpi {q) yields the free 
energy —kBT In Z for two parallel "ct" rods. Discarding 
d- independent terms, we obtain the Casimir interaction 



ifcBT^lndet [M{d, q)] 



kBTL f dq 
2^ 



In 



T{d,q) 



r(o,(?) 



'id); 



(10) 



°°dx^ 

7:— In 

27r 

4k d 
k^JI' 



1 - 



1 + d/x 



, (11) 



(12) 



It can be shown that this integral expression of the 
Casimir-like interaction is consistent with eqs. (4) and (8) 
of rcf. [S]. Our method for calculating the Casimir-like 
interaction, however, is more direct, and we shall extract 
explicit forms for the asymptotic behaviors and extend it 
to other types of rods. 

Two asymptotic Casimir interaction regimes occur 
(fig. [2]). First, in the "soft" limit where d tends to in- 
finity {k± 0), the integral in eq. ([TT|) can be performed 
analytically: 



IkBT {Ki^/KfL 



(13) 



1287r d3 

is quite small, how- 



This interaction, which scales as d^^ , 
ever, since it is proportional to d~'^ . 

Conversely, in the "hard" limit where d tends to zero 



(k± - 

^ct|ct. 



00), (jf^^"^^ reaches a constant limit, denoted by 



■ct I ct 



(d) 



^ct|ct 



'dx , 
— In 
27r 



-0.46237. 
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(14) 



(15) 



This "hard" limit corresponds to rods imposing hyy = 0. 
Note that for = [3k± — >■ 00, the q'^ factors in T{d, q) 
and r(0, g) cancel each other out in eq. pU)) . Hence, im- 
posing hyy = on the rod is mathematically equivalent 
to setting h = for x = xo and y S [0,L] (by means of 
external forces and torques). This indicates that the rel- 
ative tilt of the rods is effectively frozen in this regime. 
Physically, this comes from the fact that we are in the 
d <C i regime: a relative tilt of these long rods would 
cost a lot of energy. More formally, as we use periodic 
boundary conditions with period L in the y direction in 
our calculation, which amounts to neglecting the effect of 
the edges of the rods, hyy{xi,y) = for y G [0, i] entails 
that h{xi,y) is constant. Setting these constants to zero 
is legitimate since their presence simply corresponds to a 
possible change of reference plane. 

For ordinary rods with which corre- 

sponds to ip — ^j| ^ £±, it is the "hard" limit that is the 
physical one. Indeed, d ^ L gives d <C 4:I3kL/£± < 1, 
since the condition L < £±kBT/ {Ak) ~ 10~^^_l is required 
to avoid the in-plane undulation instability for small d/ L, 
as previously discussed. Hence, the "hard" interaction, 
eq. ()14p . is the physical one for hard rods at separations 
d <C i. This asymptotic regime is attained slowly, how- 
ever, so that one might rather expect to see an intermedi- 
ate effective power law (see fig.[21)- In principle, the "soft" 
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Table 1: Universal amplitudes A"^^ of the Casimir interactions 
between the various types of rods, rounded off to the third 
decimal place. The exact value of is — 7r/12. 
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ct ct / ct ct 


ct 
ct / ct 
ct 


-0.262 -0.357 -0.549 
-0.462 -0.672 
-0.924 



Fig. 2: Logarithmic plot of the Casimir interaction (|10p be- 
tween two "ct" rods of length L and separation d, calcu- 
lated analytically under the assumption d/L <^ 1, for L = 
= lO^^^x and k = 25 ksT (yielding d = O.ld/L). 



procedure is in agreement with zeta-function regulariza- 
tion. We thus obtain the Casimir interaction between a 
rod of type a and a rod of type (3: 



In the right shaded area the condition d/L <C 1 is violated, with A°'^^{= A^^"') given by 
and in the left shaded area the in-plane undulation instability 
is expected to occur. The dashed lines show the asymptotic 
regimes (|13|l and (|14|) . In the white region, one would measure 
PJ^^"^ Ki -0.2 (L/d)^-^. As discussed in the text, the regmie 
d ^ 1 could be compatible with d/L ^1 for anisotropic rods 
with a very soft out-of-plane bending rigidity. 



(19) 



(20) 



limit can be reached by strongly anisotropic rods with a 
large k\\ and a small kj^, such that l\_k^T / (4k) ^ d ^ L. 
Such rods should rather be called ribbons. 

Let us now calculate the Casimir interactions between 
the other pairs of rods in the d L regime. Since the 
"hard" limit is more realistic, we shall consider the situa- 
tion where the conditions hyy = 0, h^y = and h^.^^ = 
are strictly imposed. Let us define a vectorial differential 
operator D" for each rod type a £ {ct, ct, ct, ct}: 



jctlct(^) 




jctlct(^) 


= (l + 2x + 2a;2)e-2=^. 




= (l + 2a;)e-2^ 


jctlct(^) 


= 2(l + 2x2)e-2-_e-4- 


jct|5t(^) 


= 4xe-^^ + e-^\ 




^ 2e-2^-e-4^. 



{dl d^dy, 



(16) 



The constraints imposed by a rod of type a can be written 
as D"/i(x,y) = 0, for x e {0, d} and y e [0,i]. Following 
the steps detailed above in the case of two "ct" rods, we 
obtain: 



J-^l'^(d) = ifesT^lndet 



Af"l^(d, q) 



where Af"''^ is the Hermitian block matrix: 



G"l"(0,g) 


G"l^(-d, q) 


G^l"(d,g) 


G^l^'(0,g) 



(17) 



(18) 



The function f"^^ depends on the types of rods at stake: 

(21) 
(22) 
(23) 
(24) 
(25) 
(26) 

We find that the results involving "ct" rods are strictly 
identical to those involving "ct" rods: /'^'l'^ = jct|^ ^^j. 

The numerical values of the A"'^ coefficients are given 
in table [H They are all negative, which means that all the 
rods attract each other. Casimir interactions are ordinar- 
ily attractive for like boundary conditions and repulsive 
for unlike boundary conditions [3]. Here, however, the sit- 
uation is not so simple, since all the rods share the hyy = 
condition, while they differ on the other ones. We obtain 
increasing interactions upon going from "ct" to "ct/ct" 
then to "ct". Besides, as mentioned previously, "ct" and 
"ct" rods behave equivalently. However, this equivalence 
holds only in the "hard" limit. Indeed, if we put a finite 
rigidity k_l on hyy, as in eq. ([5]), the "ct" rod vanishes al- 
together in the limit — > while the "ct" rod becomes a 
different object: a flexible "ribbon" setting only h^-j^ = 0. 
We have calculated the Casimir interaction between this 



in which (G^lf^) 



G. Note that when 



type of ribbon and a "ct" rod imposing hy 



0. We 



con- find A^'M'^'^ ~0.012>0,i.e., these antagonistic boundary 



tains (e.g., for two "ct" rods), the corresponding ele- 
ment of G"l"(a;,(7) includes a term proportional to 8{x). 
This term is discarded in a; = 0, as in Ref. [17]. It is 
shown further on, in an example, that this regularization 



conditions (hyy = vs. hxx = 0) indeed yield a repulsion. 

Exact calculations of the fluctuation-induced interac- 
tions between rods at short separations d <^ L were made 
possible by neglecting boundary effects, which results in 
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Fig. 3: Logarithmic plot of tiie universal Casimir interaction 
between two rigid "ct" rods of lengtli L and separation d, calcu- 
lated numerically in all regimes of d/ L. Circles (resp. squares) 
corresponds to a length L of the rods ~ 75 (resp. 50) times 
the microscopic cutoff. The straight line shows the analyti- 
cal asymptotic behavior T^^^"^ / [k^T) = A"^^"^ d/L, valid for 
d/L < 1. 



translational invariance. Formally, this was implemented 
through periodic boundary conditions. At larger separa- 
tions, the effect of the edges of the rods cannot be ne- 
glected anymore. However, the Casimir interaction en- 
ergy can be computed numerically for all d/L. Numerical 
schemes based on the Li and Kardar method have been 
developed in Refs. [18l[T9]. The idea is to discretize the 
functional integral over the auxiliary fields ipi that enforce 
the constraints. Using a discretization scheme that will be 
discussed elsewhere [ID], we have calculated the Casmir 
interaction between two rigid "ct" rods at any separation 
(fig- El) • Although the microscopic cutoff wavelength is ex- 
plicitly involved in our calculations, we obtain universal 
results that only depend on d/L. At large separations, we 
recover the (L/d)'^ power-law obtained by Golestanian et 
al. [7]. At small separations, we recover our result (HH). 
As can be seen on fig. [31 the Casimir interaction potential 
reaches —k^T for d/L ~ 0.2 (while using the asymptotic 
form would yield the overestimated value d/L ~ 2). Thus, 
two rods approaching each other at separations less than 
half their length will effectively attract each other due to 
the Casimir effect. 

In order to check the regularization procedure men- 
tioned after Eq. ([TS]). we have recalculated the interac- 
tion between two rigid rods of the "ct" type by using 
zeta- function regularization [3T], which is rather straight- 
forward for such rods. To obtain the Casimir en- 
ergy J-^^'^^{d), we calculate the partition function Z = 
J Vh exp[-/3?^] of the membrane patch V = [0, d] x [0, L] 
situated between the rods (the free energy of the rest 
of the membrane is independent of d). The functional 
integral in Z runs over the functions h{r) that comply 
with the boundary conditions imposed by the rods, i.e., 
D'=*/i(a;,?/) = for a; e {0,d}, and y e [0,i]. Let us as- 



sume, in addition, that h{x, y) = on the rods (it will not 
change the result since we have seen previously that the 
relative tilt of the rods is effectively frozen). Finally, let 
us take periodic boundary conditions with period L in the 
y direction. Then, H = dr hOh, where O = is 

the biharmonic operator on V with the above-mentioned 
boundary conditions. Since O is positive definite Hermi- 
tian, we can write Z oc 1/ y/dct{0). Thus, discarding d- 



independent terms, we obtain J-^ 
The eigenvalues of O arc 



\d) = ifcBTTrln(0). 



n 

If 



~1X 



n e N*, me Z, 



(27) 



and the associated eigenfunctions are siii{mrx/d) x 
exp{im2TTy / L) . Hence, we obtain 



71 = 1 



dm \n{Xnrn) . (28) 



While Eq. (|28p is divergent, its finite Casimir part can 
be extracted by calculating Co{s) = S -oo'^'^ Kirm 
which is finite for s > |;, and yields by analytical contin- 
uation: 



•ct I ct 



(d) = -ifcBTlimC^(s) 

Z s— >0 



TT L 

-'-^12 d 



(29) 



This result is in agreement with the value A'^*!'^' = — 7r/12 
obtained above with the regularization of Li and Kar- 
dar [17| . This agreement also implies that a half-plane 
imposing the same boundary conditions as a rod will yield 
the same Casimir interaction (in the regime L ^ d). Phys- 
ically, this equivalence is due to the fact that the Casimir 
interaction arises from the confinement of the fluctuation 
modes between tire rods (or the half-planes), i.e., in the 
patch v. 

Finally, we have checked that the undulation instability 
discussed in Ref. |8] can be recovered from the Casimir 
interaction calculated in this paper. Let us consider the 
fundamental deformation mode: 



Sx^u{y)^Usm{ny/L'), ye[0,L'] 



(30) 



where 6x is the in-plane displacement of each rod in the 
direction of the other rod. Inextensibility implies L' — 
L - tt'^U'^ /{AL) + 0{U^). In the limit < L, wc may use 
the proximity-force (or Derjaguin) approximation j22j to 
calculate the total free energy of the two rods: 



cf 



L' 



dy 



k^T Jq \d 
d nd 



2u 



nU 
'2d 



d 



(31) 



Minimizing J-def with respect to U , we find that the equi- 
librium deformation amplitude [/eq reaches d/2, i.e., that 
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the two rods touch each other when d < d' where 



L 



-4(2 + 7r)A 



a\l3 



1/3 



1/3 



(32) 



We thus obtain an instabihty threshold very similar to that 
of eq. ([T]). For instance, for two "ct" or "ct" rods, with 
^«l/3 ^ -0.46, the prefactor of {L/ip)^/^ is 0.31, while in 
eq. © it is 0.27. 

In this paper, we have assumed that the effect of the 
membrane tension a was negligible with respect to that 
of its bending rigidity k. This hypothesis is valid in usual 
flaccid membranes (cr « lO^^N/m) for rods shorter than 
(K/ay/-^ « 1 ^m. However, a can be increased up to about 
10~^N/m by applying an external tension to the mem- 
brane: in this case, the effect of tension will become impor- 
tant. The methods presented here in the bending- rigidity- 
dominated regime can also be applied to the tension- 
dominated limit, where L and d are much larger than 
{n/aY^^. For instance, the interaction between two "ct" 
rods imposing hyy = is then still given by Eq. [T21 but 
with = -7r/24. This result can be obtained both 

through the method based on Ref. [17] and through zeta- 
function regularization. This tension-dominated interac- 
tion is relevant in liquid- vapor interfaces [23]. A general 
calculation of the Casimir-like interaction between two 
rods, including both k and a, is beyond the scope of this 
paper. This problem has been addressed very recently in 
the simpler case of two small circular inclusions [24] . 

In summary, we have studied the Casimir interaction 
between two parallel rods adsorbed on a fluid membrane 
characterized by its bending energy and subject to thermal 
fluctuations. We have considered four types of boundary 
conditions, depending on whether the membrane can twist 
along the rod and/or curve across it. Whatever their type, 
rigid rods of length L at short separations d <^ L, experi- 
ence an attractive Casimir interaction scaling as kBTL/d. 
In this regime, two of the four types of rods are equivalent, 
which yields six universal Casimir amplitudes. When tak- 
ing into account the finite persistence length £p of the rod, 
we find, in agreement with Ref. [8] , that at separations less 
than dc ~ L{L/(.pY/'^ two rods with fixed ends will bend 
toward each other and finally come into contact because 
of the Casimir interaction. Besides, we have shown that 
rods with a very soft out-of-plane bending rigidity, but 
a large in-plane bending rigidity, still experience a power- 
law Casimir interaction (albeit small), which scales as d~^ . 
For a large, but not infinite out-of-plane bending rigidity, 
an effective power law in d"" with a close to 1 would 
be observed. This nuances the conclusion of Ref. [5] ac- 
cording to which a finite out-of-plane rigidity, no matter 
how large, should destroy the long-range Casimir inter- 
action. We have also calculated numerically the Casimir 
interaction between two rods at any separations, finding 
at d/L « 0.5 a crossover between the short separation 
c?~^ behavior and the large separation d"^ behavior. We 
have shown, in the case of the simplest boundary condi- 



tions, that at separations d/L < 0.2 the Casimir attrac- 
tion reaches k^T. In conclusion, due to the very long 
range of the Casimir interaction, long rods adsorbed on 
an effectively tensionlcss membrane will actually attract 
each other — by the Casimir effect — when their separation 
becomes significantly smaller than their length. 
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